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DETERMINATION OF STRESSES IN ELASTIC SOLIDS USING 


THREE STRESS FUNCTIONS AND THREE EQUATIONS 
By Robert E. Reed, Jr. 

Ames Research Center 

SUMMARY 


The classical equations of elasticity in terms of the stress components 
are reduced to a set of three equations in terms of the three Maxwell stress 
functions and arbitrary functions of integration. The sufficiency of these 
equations is due to the interrelationships between the six compatibility equa- 
tions. The difficulty in choosing the arbitrary functions to fit a particular 
problem accents the advantages of using the inverse method to find solutions. 

A family of exact solutions is presented for stresses in a rectangular solid 
with certain applied stresses on two opposite faces, workless boundary condi- 
tions on two opposite faces, and the remaining two faces free of applied 
stress . 


INTRODUCTION 


The differential equations for determining the stresses and displacements 
in a three-dimensional elastic solid are well known (e.g., refs. 1 - 3 ), but 
relatively few exact solutions have been found. Primarily, two methods can be 
used to find these solutions . One method is to choose displacements that 
satisfy the equilibrium equations and boundary conditions and the other method 
is to choose stress components that satisfy the equilibrium equations, bound- 
ary conditions, and compatibility equations. The choice between these two 
methods will depend on the quantities desired and the type of boundary condi- 
tions present. However, this report will deal only with the method of choos- 
ing the stress components; and then it attempts to derive a simplified set of 
equations from which exact solutions can be obtained. Basically, the complete 
system of equations for the stress components consists of 12 unknowns and 
15 equations. The unknowns are the 6 stress components and the 6 strain com- 
ponents. The equations consist of 3 equilibrium, 6 stress-strain, and the 
6 compatibility equations in terms of the 6 strain components which are neces- 
sary and sufficient to ensure the continuity of the displacements (for simply 
connected regions) of any solid undergoing deformation. This system can be 
reduced to 6 equations and 3 unknowns by defining 3 stress functions that 
identically satisfy the equilibrium equations and then, with the stress-strain 
relations, by introducing these stress functions into the compatibility equa- 
tions. However, it is pointed out in some references (e.g., refs. 1-3) that 
the 6 compatibility equations are interrelated and that, at most, 3 of the 
6 are unrelated equations. That is, certain groups containing 3 compat- 
ibility equations will have no equation related to the other 2 , but any group 



of k or more equations will have 1 or more equations related to the other 
equations in the group (see appendix A for a discussion of the compatibility 
equations). These interrelationships, and the assumption that the material is 
homogeneous, isotropic, and follows Hooke* s law, make it possible to derive a 
set of 3 equations in terms of 3 unknown stress functions and some arbitrary 
functions of integration. These equations, along with the boundary condi- 
tions, are necessary and sufficient for determining exact solutions, but the 
arbitrary functions of integration must be chosen to suit each particular 
problem. Exact solutions of these equations are often difficult where a par- 
ticular problem is specified. However, use of the "inverse method," where the 
arbitrary functions of integration are chosen such that solutions can be found; 
will produce exact solutions to problems which may be of immediate interest in 
terms of specific problems or can be used to check the validity of approximate 
solutions which are based on formulations that can be used to solve a wider 
range of problems . 

In the following, the set of three equations is derived and an example 
solution is presented and discussed. This solution represents a rectangular 
solid with certain combinations of shear and normal stress applied to opposite 
faces while two other faces have workless boundary conditions and the 
remaining opposite faces are free of applied stress. 

SYMBOLS 


Ai 

a 

B± 

b 

E 

F x^ F y^ F z 

fjgjkA 
F,G,H J 


i 

Lo 


constants of integration 
dimension in x direction 
constants of integration 
dimension in y direction 
elastic modulus 
resultant applied forces 

functions of integration 

a 

b 


P,Q,R 

P i'W U i 

u,v,w 


resultant applied moments 
functions of body forces 
functions of | 

displacements in x,y,z directions 
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I 



UjVjW 

X, Y,Z 


P,7,A 

7xy>7 xz >7y Z 

€ x> € y> e z 

0 U@2 >&3 
V 

S*ib£ 

a x'V ff z 

T xy' T xz > T yz 


dimensionless displacements; X , X and — 

w* D Si 

"body forces 

rectangular coordinates 
constant parameters. A 2 = 7 2 - P 2 L 0 
shear strains 
normal strains 

functions of >?>Q> and R 

Poisson* s ratio 

dimensionless coordinates; — , ^ , and — 

ah a 

normal stress components 

shear stress components 
Maxwell stress functions 


Subscripts 

I imaginary part of a complex quantity 

R real part of a complex quantity 

HOMOGENEOUS, ISOTROPIC, LINEARLY ELASTIC MATERIALS 


The equilibrium equations including the body forces are 


Scr St 


Sx 


+ 




St 


xz 


dy 


Sz 


+ X = 0 


do y ^ ^ T xy + ^ T yz + 


dy 


Sx 


Sz 


Y = 0 > 





+ Z = 0 

J 


( 1 ) 
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Hooke 1 s law is 


G x g v ( a y + o’z ) ] 

€ y = jg [ °y ~ V ( °X + °Z ^ 

€ z - — [cr z - v(a x + Oy) ] 

_ 2(1 + v) _ 2(1 + v) 

^xy ^ Tx y Jl ^xz T xz ' 




\ ( 2 ) 
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E 


yz 
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and the compatibility equations are 
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^ 7 y Z 

d7 xz 

+ s = 

dy 2 

dx 2 

dx by 

L— 

by dz 
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^ 7 xy\ _ 
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dy 
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(3) 


The equilibrium equations can be satisfied by the Maxwell stress func- 
tions modified to account for the body forces* Denoting d 2 (»«)/dx 2 "by 
( # *),xx> e ^ c • y one obtains 


°x = jZZ + * 3 ,yy " J* X dx ' T xy = "*3, 


xy 


T xz ~ 


2 ,XZ 


0 = + \|f - f Y dy y 

y ljyZZ 3^XX J J 

= >1 f + \1/ - P Z dz , T Trr7 = - 1 |r 

z v i,yy s^xx J ' yz i,yz 
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Equations (k) can be substituted into equations (2) to give the strains in 
terms of the stress functions. Equations (2) can then be substituted into the 
compatibility equations to give six equations in terms of the three unknowns 
\|f \|r , and \|r 3 . For convenience, the following definitions are made: 


01 = w 2 ^ - V% - v 2 t 3 + (t ipcx + t 2jyy + % )ZZ ) + P 
e 2 = -v 2 t x + w% - v z % + (i 1 ^ xx + * a ,yy + t 3 ^ zz ) + Q 


0 3 = 

-v 2 ^ - v% + w% + (t ipDC + + 

^3 ,ZZ^ 

where 

P = J* X dx - v(j Y dy + J 

N 

' Z dz 
/ 

) 


Q = J Y dy - v(j X dx + J* 

\ 

' Z dz 
/ 

) 


R = J Z dz - v(j X dx + J 

\ 

' Y dy 

/ 

) 

and 

n2 >2 

V 2 (..) - - (..) + (••) 

dx 2 Sy 2 

+ 

dz 2 

(••) 

Equations (3) then become, respectively (see appendix B for 
derivation) , 


> 




0 i,yy + 0 2 ,xx = 0 




0 1,ZZ + ®3,XX = ® 

> 



R 


( 5 ) 


0 2,zz + 0 3,yy 0 


(6) 


and 


e i,yz 0 


2 ,XZ 


= 0 > 


0 3,xy = 0 


( 7 ) 


l 



Equations (7) can*be integrated to find 9 i, 9 S , and 0 3 , but these solutions 
must be restricted to satisfy equations (6). Therefore, equations (7) yield 


01 = f i(x,y) + gi(x,z) 

0 2 = f 2 (x,y) + h. 2 (y,z) 
63 = g3(x,z) + h 3 (y,z) 


(8) 


and equations (6) require that 


fi.yy + f z,xx = 0 


Sl,ZZ + £>3, XX ^ 


> 


h 2 ,zz + h 3,yy 0 


(9) 


Equations ( 9 ) can be integrated to give 


f i( x ,y) = -JJf 2 ( x ^y),xx d y d y + y f 3 ( x ) + F ( x ) 

g 3 (x,z) = -J‘Jg 1 (x,z)^ zz dx dx + xg 4 (z) + H(z) ) 

h 2 (y, z) = qjh 3 (y,z) dz dz + zh 4 (y) + G(y)^, 


( 10 ) 


Substituting equations (10) into equations (8) gives 

= -JJ’f 2 ( x ,y), xx dy dy + yf 3 (x) + g x (x,z) + F(x) 

0 2 = _ JJ h 30 ^z),yy <3- z <3- z + z h 4 (y) + f 2 (x,y) + G(y) \ 

6 3 = ~H s ^ x ’ z Kzz dx dx + ^4(2) + h 3 (y^ z ) + H(z) 

j 


(ID 


Equations ( 5 ) can be used to write equations (ll) in terms of the \|r values. 
This gives the three equations which, along with the boundary conditions, 
form a necessary and sufficient set for determining elasticity solutions. 
These equations are 
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W 2 ^ - V% - V% + (t 1;XX + * 2 ,yy + *3,BZ> + P 


"\ 


= -J'f f 2 ( x ^ ) ,xx d y d y + yf 3 ( x ) + 6iU> z ) + F U) 


-v 2 ^ x + - v% + (* 1}XX + % >vv + i 3)ZZ ) + 


Q 


= -J’J’h 3 (y,z)^ yy dz dz + zh^Cy) + f 2 (x,y) + G(y) 


> ( 12 ) 


- V 2 ^ 2 + W 2 i 3 + (t^XX + * 2 ,yy + f^ zz ) + R 

= -JJsi(x,z)^ zz dx dx + Xg 4 (z) + h 3 (y,z) + H(z) 


J 


where the right-hand sides are arbitrary functions of integration. Note that 
there are no arbitrary functions of the form f(x,y,z). 

Particular solutions can be found for F(x) , G(y) , H(z) , and yf 3 (x) , 
zh^Cy), xg 4 (z); they have the form: 

iip = 1 ^ - iitF(x) + yf 3 (x)]dx dx 


= 2 . +V + zh 4(y)3 d y d y 


t 3p = 3^7 JJ( H ( Z ) + xs 4 ( z )3 dz dz 

The stresses resulting from these solutions are zero so these functions can 
be disregarded. However, the solutions for f 2 (x,y) , g x (x,z), and h 3 (y,z) 
will yield nonzero stresses which may be necessary to satisfy the boundary 
conditions of a specified problem. 


COMMENTS ON POSSIBLE SOLUTIONS 


General Comments 


The relevant arbitrary functions (f 2 , g x , h 3 ) in equations 
expanded into appropriate double Fourier series, and particular 


(12) can be 
solutions can 
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be found in terms of the constants of the series. These constants would be 
determined to fit a particular problem. However, these functions and the 
body forces will be taken to be zero and only the homogeneous equations will 
be considered. These are 


W 2 t i - V% - V% + (* XtX Z + + % }ZZ ) = 0 

-V 2 *! + VV% - V% + (tl,xx + * 2 ,yy + * 3 _ >zz ) = 0 > 


(13) 


Subtracting equations (13) 


in pairs gives 


v 2 i 1 - v 2 f 2 
= V 2 ^ 3 , 


(14) 


Therefore, any solution of equations ( 13 ) will satisfy equations (l4) . 
This gives some indication of the type of solutions one might choose. Another 
possibility is to let 


= ^2 = ^3 = * 


Equations (13) then reduce to a single Laplace equation: 


V 2 ^ =0 (15) 

From equations (4), it is seen that solutions to equation (15) represent the 
case when the first stress invariant (a x + + a z ) is zero. One obvious 

form of solution to equation (15) is ^ = A cos ax cos Pye^ z , which is 
similar to a solution discussed in reference 4. 


The solution to equations ( 13 ) discussed here gives the stress distribution 
in a rectangular solid having two opposite sides free of applied stress, two 
opposite sides constrained by workless boundary conditions, and the remaining 
two sides subjected to some combination of applied stress. 


Solution to equations ( 13 ) The following dimensionless coordinates will 
be used: 


t = £ 

* a 9 


T = ? 


t = 


a 


The term V 2 (") is then given by 

a 2 V 2 ( * • ) = (■•),!! +Lo(-'), T jt 1 + (••),{;£ 
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The solution is chosen in the form 


\ 

*3. = f(6)y(ri)z(c) 

= g(|)Y(ri)Z(5) > 
= h(|)Y(ri)Z(0 

where it is further assumed that 

trn = -P 2 ^ > Z ,£C " 


(16) 


(17) 


Letting A 2 = 7 2 - P 2 L 2 > one can separate the variables in equations (13) by 
substituting equations (l6) . The resulting ordinary differential equations 
are 

(1 + v)f" + vA 2 f - (g" + 7 2 g) - (h" - ^lJi) = 0 ( 18 ) 

-A 2 f + [vg" - (1 + v ) P a L 2 g + V7 2 g] - (h" - P 2 L 2 h) = 0 (19) 

-A 2 f - (g" + 7 2 g) + £v(h" - P 2 L a h) + (1 + v) 7 2 h] = 0 (20) 

where f M = (d 2 f/d£ 2 ). Subtracting equations ( 19 ) and (20 ) gives 

g" + A 2 g = h ,f + A 2 h 


so 


g = h + Ai sin A£ + A 2 cos A£ 


( 21 ) 


In the following, Ai and A 2 are assumed to he zero (which would he required 
later to satisfy the boundary conditions) . Substituting equation (21 ) into 
equation ( 19 ) or (20 ) gives f in terms of h. This relation and equa- 
tion (2l) can be substituted into equation (l8) to give the following equation 
for h(§); 

h ,m + 2A 2 h ,f + A 4 h = 0 (22) 

Equation (22) has the solution 

h(£) = B x e lA& + B 2 e~ lA ^ + B 3 iA£e lA ^ + B 4 iA£e~ lA ^ (23) 

The functions f(|) and g(^) can now be found and the stresses can be 
determined from equations (4) . 
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The desired boundary conditions of two opposite surfaces that are free 
of applied stress are given by: 

|=0 


°x T xy “ T xz - o 


1=1 


a x T xy _ T xz “ 0 


Since and t xz are both zero when h’(|) = 0, these six conditions give 

the following four homogeneous equations: 


h(o) = 0 
h'(0) = 0 
h(l) = 0 
h'(l) = 0 j 


(24) 


For the B values to be nonzero, A must have such a value that the determi- 
nant of the coefficients in equations (24) is zero. The B values can then 
be found in terms of one arbitrary constant from equations (24) and 
substituted into equation (23) to give 

h(| ) = a 2 B Q A[(l - f-)sin A sin A£ - A£ sin A(l - |)] (25) 

In order for the determinant to vanish, A must satisfy one of the two 
equations 


sin A = ±A (26) 

The nonzero roots of equations (2 6) are complex and the lowest five roots of 
each of equations (26) are given in table I (see, also, ref. 5)- 


TABLE I.- ROOTS OF sin A = ±A 


sin A = -A 

sin A = A 

a r 

A I 

a R 

A I 

4.21239 

10.71253 

17-07336 

23-39835 

29.70811 

2 .25072 
3 .10314 
3 .55108 

3 .85880 
4.09370 

7.49767 

13 -89995 
20.23851 
26.55454 
32 .85974 

2.76867 

3-35220 

3 .71676 
3-98314 
4.19325 
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If equation (25) is used, but not equations (26), which would require the 
± sign, the stresses obtained will be complex quantities given by 

cr* = A s B 0 [(l - | ) sin A sin A| - A| sin A(l - |)]Y(ti)Z(^) 

ct* = B q |aP 2 Lq[( 1 - | ) sin A sin A| - A£ sin A(l - £ ) ] 

+2(0 2 L^ - V7 2 )[sin A cos A| - A cos A(l - | ) ]j-Y(r))Z(£) 

ct* = B 0 |^7 2 A[-(l - |)sin A sin A| + A| sin A(l - |)] 

+2(7 2 - v0 2 L 2 )[-sin A cos A| + A cos A(l - | ) ]j-Y(t])Z(£) 

* / (27) 
= AL Q B 0 [sin A sin A| + A sin A(l - |) - A(l - |)sin A cos A| 

-A 2 1 cos A(1 - | ) ]Y' (t])Z(£) 

t* = AB 0 [sin A sin A| + A sin A(l - |) - A(l - |)sin A cos A| 

-A 2 | cos A(1 - |)]Y(ti)Z*(5) 

t* z = L 0 B 0 -j^-A(l - |)sin A sin A| + A 2 | sin A(l - |) 

-2(1 - v)[sin A cos A| - A cos A(l - | ) ]j"Y' ( rj ) Z * (£) j 

Both the real and imaginary parts of equations (27) are solutions to the equi- 
librium equations, boundary conditions, etc., so the actual stresses can be 
taken as the real parts of equations (27)* When equations (26) and the 
notat ion 


7 = 7 R + i7j 
B q = % + iBj 

Z(5) = e 7 ^ = e 7R ^(cos 7 j(; + i sin 7i0 
are used, the real stresses are 
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o x = -Y(.T))e 7R ^B R [P 1 (|)cos 7 X ^ - P 2 (|)sln 7jO 
-Bj[P 2 (^)cos 7 j^ +P 1 (|)sin 

Oy = -Y(Ti)e 7R ^B R [P 3 (|)cos 7jC - R 4(t)sin 7j<L] 

-Bj[P 4 (5)cos 7j^ +P 3 (|)sin 7iU} 

a z = Y(n)e 7R ^B R [P 5 (|)cos y-£ - P 6 (|)sin y-£] 

-Bj[P 6 (| )cos 7 t C +P 5 (|)sin 7 t^]| 1 

J > (28) 

T xy = L 0 Y , (il)e 7R ^{B R [Q 1 U)cos 7j£ - Q 2 (|)sin 7jU 

-Bi[Q 2 (a)cos 7jC +Q 1 (|)sin 7 j£]} 


T xz = Y(ri)e 7R ^B R [Q 3 (i)cos 7 ^ - Q4U )sin 7jU 
-B I [Q 4 (a)cos 7j^ + Q s (i ) sin 7 j5]} 

T yz = L 0 Y'(Ti)e 7R ^B R [Q 5 (i)cos 7j i ~ Qs(l)sin 7^] 
-BjtQeUJcos 7 i£ + Q 5 (i )sin 7^]} j 


where Pi(|), . . . , Q 6 (|) (defined under Definition of Functions) contain 
the ± sign of equations (26) . 


The displacements can he found by integrating the strain-displacement 
equations: 


8u 8v 



8w 


The arbitrary functions resulting from the integration are taken to be zero to 
satisfy the shear strain-displacement equations. The dimensionless 
displacements are given by 
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U = - Y(r 1 )e 7 R ^B R [U 1 (|)cos 7 Z t - U 2 (|)sin 7 ^] 

-Bi[U 2 (|)cos 7jC + Ui(g ) sin 7 i£]} 

v = - 1 * - P 2 L^Jy dr)je 7:R HB R [U3(5)cos 7-^ - U 4 (| )sin 7-,-^] 
-B i [U 4 (|)cos 7 Z C + U 3 (|)sin 7^]} 

w = — jp Y(T 1 )e 7R ^B R [Q 5 U)cos 7 ji - Qs(i)si n 7 Z U 
-BjtQsCl )cos 7 j £ + Qs(| )sin 7 iU} 


(Ui(|), . . . , U4(|) are defined in the next section). 

The functions Pi(g), . . . , U 4 ( ^ ) depend on which of equations ( 26 ) is 
chosen and, depending on the choice, the resulting stresses and displacements 
will have different symmetry characteristics. For 

sin A = A : 

a x’ a -y’ a z> r y Z >v>w are odd about £ = g, that is, cr x (£,r),£) = -cr x (l - £,tj,{;) 
T xz’ T xy’^ are even about £ = |, that is, t xz (|,t],^) = t X2 (i - 


sin A = -A: 


a x> a y> a z> r yz> v ’™ are even - about £ = ^ 
T xz * T xy 9 u are °dd about £ = i 


Since A is determined from one of equations (26) and the relation 
A 2 = 7 2 _ p 2 L^ must be satisfied, only one of the three parameters 7, 0, 
and A remains arbitrary. In the following, the parameter 0 and Y(r|) are 
chosen to satisfy boundary conditions at T) = 0, 1. Two obvious choices are 

Y = cos 0tj 0 = njt , n = 0,1,2, . . . 

V = T X y ~ T yz ’ 0 a *^ T 0,1 


and 
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Y = sin |3 t] 0 = nn , n = 1,2, . . . 

Cy = u = w = 0 at T] = 0,1 

Since each, of equations (26) contains infinitely many roots and n can 
have any integer value, the stresses and displacements can be summed to give, 
for example, 

° x = - ^ cos M e ^ {BRmntPim(l)cos 7^5 - P 2m (|)sin 7 Trrin U 

m=i n=o 

-B Tmn ( | ) COS 7 Tmn £ “ ) s ^- rl 7 Tmn C Jj* 

Two sets of arbitrary constants are available to satisfy boundary conditions 
on the surface £ = constant , so one could specify two stress components on 
£ = £ 0 and would have to accept whatever the solution gave for the other 
stress applied to the £ = constant surface. The functions cos 0 n t] are 
orthogonal so a function can be easily represented in the T} direction; but 
the functions P m (^) , etc. , are not orthogonal so it could be difficult to 
approximate a specified stress in the | direction. To avoid this difficulty 
in the present report, the numerical work is confined to single-term solutions 
given by equations (28). 

Definition of functions .- The functions that' comprise equations (28) and 
(29) are defined in this section. First, the constant parameters are 

7 2 = \ \J( p 2 Lq + ^i) 2 + *2 - (p 2l ? + 


A r Aj 


Ai = a; 


A 2 2A r Aj 


The functions that depend on the choice of equations (26) are 


P R (|) = T (l - |)sin A R | cosh Aj£ + | sin A R (l - |)cosh Aj(l - |) 
Pj(|) = + (l - |)cos A r | sinh Aj£ + £ cos A R (l - £)sinh Aj(l - £) 


Q R (| ) = +cos A R £ cosh Aj£ + cos A R (l - £)cosh Aj(l - £) 
Ql(i) = ± sin A r | sinh A-j-£ - sin A R (l - £)sinh Aj(l - |) 
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T R (|) = ±sin Ar£ cosh. Aj£ + sin A R (l - |)cosh Aj(l - |) 
T j(^) = ±cos A R | sinh Aj^ + cos A R (l - g)sinh Aj(l - | ) 


U R ( | ) = +(l - 5 )cos A R | cosh Aj| - | cos Ar (1 - |)cosh Aj(l - 5 ) 

Uj(|) = ±(l - |)sin A R £ sinh Aj^ + i sin Ar( 1 - |)sinh Aj(l - |) 

where the sign corresponding to that of equations ( 26 ) must he used (i.e. 

±A = sin A) . The functions appearing in equations ( 28 ) and ( 29 ) are 

Pi(l) = (Ai - A 2 )P R U) - SAiAaPjU) 

Pa(i) = 2 AiA 2 Pr(|) + (A* - A|)Px(|) 

P 3 (l) = P 2 Lq[A!P r (|) - A 2 P i (|)] + &[a r [(1 - v)p 2 L l - vA a ] + vAjA 2 }q r (|) 
-2 -[a i [(1 - v ) (3 2 Lq - vAj - vA r A 2 }q i (|) 

P 4 (|) = P 2 L 2 [A 2 P R (|) + A 1 P 1 U)] + 2-[a i [(1 - v)P 2 L 2 - vA a ] - vA R A 2 j-Q R (| ) 
+2 -{Ar[( 1 - v)p 2 L 2 - vAj + vA I A 2 j-Q I (|) 

P 5 (i) = (At + Ao.L 2 p 2 - A 2 )Pr(|) - A 2 (P 2 L 2 + 2Ai)P I (|) 

+2{a r [A! + (1 - v ) |3 2 L 2 ] - A I A 2 j-Q R (| ) 

-2^Aj[Ai + (l - v ) P 2 L 2 ] + A r A 2 }qi(|) 

P 6 (|) = A 2 ( 0 2 L 2 + 2Ai)Pr(|) + (A 2 + A X P 2 L 2 - a|)Pj(|) 

+2"|^Aj[Ai + (l - v)p 2 L Q ] + A R A 2 ]-Q R (|) 

+2{a r [A x + (1 - v)p 2 L 2 ] - A i A 2 }q i (^) 

Qi(t) = AiT r (^) - A 2 Tj(| ) + (A r Ai - AjA 2 )U r (|) - (AiAq + A r A 2 )Uj(|) 

Q 2 (i) = A 2 T r (|) + AiTj(l) + (AjAi + A r A 2 )U r (|) + (AiA r - A 2 Ax)Uq(|) 

^ 3 ( 1 ) = (7rAi - 7jA 2)T r (|) - (/jAi + 7 rA2)Tj(^) 

+[Ai(7rAr - 7lAj) - A 2 (Ar7j + Ax/r)]Ur(|) 

-[Ai(7rAx + AxAr) + A 2 (Ar7r - Ax7i)]UxU) 



Q 4 (|) = (7jAi + 7 R A2)T r (|) + (7 R Ai - 7 jA2)Tj(|) 

+[Ai(7jA R + 7 r Aj) + A2(7rA r - 7jAj)]U^(|) 

+[Ai(7rA r - 7jAj) - A 2 (7jA r + 7 R Ai)]Uj(i) 

Qs(l) = (A R Aa - 7 iA 2)P r (|) - (7 R A2 + 7 jAi)Pj(|) + 2(l - v)(Aj^7 R - Aj7j)Q r (|) 
-2(l - v)(A r 7j + Aj7 R )Qi(|) 

Q s (|) = (7 r A2 + 7iAi)P r (| ) + (7 R Ai - 7jAa)Pi(l) + 2 (l - v H7jA r + 7 r Aj)Q, r (|) 

+2(1 - V )(7 R A R - 7jA j) Qj ( | ) 

Ui(|) = (1 - 2 v)[AiT r ( 6) - AaTjd)] - (AiA r - A 2 A I )U R (|) + (AiAj + ^ArJUjU) 
U 2 ( | ) = (l - 2 v)[A 2T r (|) + AiTj(|)] - (AiAj + A 2 A R )U R (|) - (AiA r - A 2 A];)Ujd) 
u 3 (|) = AiP R d) - A 2 Pid) + 2(1 - v)[A R Q R d) - AjQjd)] 

U 4 (i) = A 2 P R (i) + AiPj(i) + 2 (l - v)[AjQ R d) + A r Qj(| ) ] 

NUMERICAL EXAMPLES 


The numerical evaluation of any of the solutions given by equations ( 28 ) 
is straightforward and is easily carried out, but only the solution for the 
lowest value of A and P = 0,n is presented here. The boundary conditions 
satisfied by the solutions considered are 


S = 0 , 1 : 


T) = 0;1 : 


-"x 


= T 


xy 


= T 


XZ 


= 0 


V - Tjjy - T y Z - 0 


£ = 0 ; 


a z = cos Pt][B r P 5 (|) - BjP e (s)] 

T xz = cos Pti[B R Q 3 (|) - B jQ 4 ( | ) ] 
x yz = -PL 0 sin Pr][B R Q 5 (i) - BjQsd)] 




oo; 


°z T xz T yz 


= 0 
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Since the condition of decreasing stress with increasing z is imposed, 7 j 
is chosen to he negative. The other parameters are 

sin A = -A 

(3 = 0,jt(P = 0 corresponds to plane strain) 

Y(tj) = cos Prj 

L 0 = 0.01, 0.5. 1.0 


Combinations of these parameters are calculated for the following two cases: 


Case 1: 

Br f 0, 

Bj = 0 

( a Xl> 

CT yi> 

etc . ) 

Case 2: 

•N 

0 

II 

pq 

Bi i 0 

(°X2 3 

a ya' 

etc . ) 

For both cases, the applied stresses on the 

surface 

i = 

0 can produce resul- 


tant forces and moments that can he found from the following relations: 


Ft = 



- " x J t=0 11 an 


“s 




C =o d5 a,) 



T yzl £_ 0 d £ dr l 




d| drj 



<*zl, d £ dT l 

£=° 




where 

F |> F r|> F £ 

M|, Mjj-, 


£5 

1 ^ 

1 ^ 
In 

ah 9 

ah ah 

ik 

My M z 

ah 2 

' a 2 b ' ab 2 


and the notation is shown in 

Sketch (a) sketch (a). For the given choice of 

parameters, the resultants are 


IT 



= 0 



Fj, = -4vL 0 [(-l) n - 1][B r (7rAr - 7jAj) - Bj(7jA r + 7 r Aj)] 

M| = -4v[(-l) - l](B R A R - BjAj) 

For plane strain (n = 0), all resultant forces and moments are zero. 

The calculated stresses for cases 1 and 2 are shown in figures 1 through 
12 and several linear combinations of the two cases are shown in figures 13 (a) 
to (d) and l4. Although some discussion of the limited numerical work is 
given here, a much closer examination of the entire family of solutions 
presented here would undoubtedly yield much useful information. 


Discussion of Numerical Results 

Some general comments should first be made about the figures . All the 
stresses were calculated at the points | = 0, 0.1, 0.2, 0-3; 0.4, 0.5 so 
the graphs are accurate at those points and the stresses are either even or 
odd extensions for 0*5 < I < 1 as indicated on the graphs. Also, all the 
graphs are independent of rj since this coordinate is contained in the 
abscissas . 

For case 1, the applied stresses are shown by the £ = 0 curves in fig- 
ures 1, 3, and 5- The remaining three stress components are shown in 
figures 7, 9, and 11. 

For each stress component the same curve applies for P = 0 and (3 = jt 
when L Q = 0.01 ( T xy = T xz “ 0 as L 0 0 or P = 0) which means that the 
stresses for P = % and L Q = 0.01 can be represented by cos Prj 

where a — are the stresses for the plane strain problem given by (3 = 0. 

In fact, very little difference is apparent up to L Q ~ 0.1 (the accuracy is 
determined by the smallness of PL Q rather than just L 0 ). This is useful 
to know since the plane strain solution to this and similar problems can be 
obtained from two-dimensional theory. 

The graphs show the rapid decay of all the stress components as £ 
increases. However, o xx requires a longer distance to reduce its value to a 
negligible amount because it has its maximum value at about £ = 0.2 whereas 
the other components have their maximum values at £ = 0. At £ = 1 (a dis- 
tance equal to the thickness a), a xi is about 5 percent of its maximum 
while the other components are about 1 to 2 percent of their maximum value. 

A comparison of the curves for L Q = 1.0 and L Q = 0.01 for p = jt shows 
that the applied stresses are similar in shape but the magnitude of each com- 
ponent is larger for L 0 = 1.0. Two of the remaining stresses, a xi and 
have similar shapes but Oj ± changes both its shape and magnitude. This 
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shows that despite similar applied stress distributions, the internal stress 
distribution may vary considerably with the thickness ratio. 

The applied stresses for case 2 are shown in figures 2, 4, and 6 , and the 
remaining stresses are shown in figures 8 , 10, and 12. The main difference 
between the applied stresses of cases 1 and 2 (figs. 5 and 6 ) is the distri- 
bution of T y Z i a nd T yz 2 * '^ ie same comments as for case 1 can be made regard- 
ing the similarity between the 0 = 0 and 0 = it solutions and the decay of 
stress with increasing £ except for the behavior of compared to that 

of cr xi . In figure 8 , one sees that has its maximum value at £ = 0, 

and the maximum value is an order of magnitude larger than that of cr xi . This 
is true for the solution 0 = 0 for which Ty Z1 = T y Z 2 = 0 so the difference 
in distribution of T yz does not explain the large maximum value of cr-^. 
Linear combinations of cases 1 and 2 can be formed to give a better under- 
standing of what governs the behavior of cr x . Considering only the solutions 
for 0 = 0 , one sees that cases 1 and 2 have the resultant forces shown in 
sketch (b). The shear forces are given by the area (0 < £ < 0.5) under the 


(* 

(* 


■z 


( 


Cose I Case 2 

Sketch (b) 


curves for T XZ1 and T XZ2 for £ = 0 (figs. 3 and 4). The moments are caused 
by a z , as shown in figures 1 and 2. The relative value of the force and 
moment can be varied by superposing cases 1 and 2 ; four examples are shown in 
figures 13(a) to (h) . The values indicated in the figures for the force and 
moment are not precise since they were determined graphically. In fig- 
ures 13 (a) to (c) the force is kept constant but the moment decreases from a 
dimensionless value of about 50 in figure 13(a) to zero in figure 13(c). The 
maximum value of o x increases approximately linearly as the moment 
decreases. Figure 13 (d) shows the case for which the force is zero and the 
moment nonzero, so figures 13 (c) and (d) show the uncoupled effect of a force 
or moment on the stress components regardless of the value of L 0 . This 
points out that it is possible even in thin plates to have large transverse 
stresses near the edge of the plate, but these stresses, in agreement with 
Saint Venant f s principle, are contained within an edge region whose width 
about equals the plate thickness. This is clearly seen in figure l4 where 
a x at £ = 0.5 is plotted versus £ for each example given by figures 13 (a) 
to (d) . All the stresses have decreased to small fractions of their maximum 
values within the distance £ = 1 . 
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CONCLUDING REMARKS 


The set of six fourth-order compatibility equations in terms of the three 
Maxwell stress functions has been -written in the form of three second-order 
equations with arbitrary functions of integration to simplify the task of 
finding exact solutions by the inverse method. The solutions studied to date 
are for body forces and arbitrary functions of integration being zero. The 
most interesting one found is presented in the form of an infinite family of 
solutions which gives the stresses in a rectangular solid subjected to applied 
stresses on two opposite surfaces with two other opposite faces constrained 
by workless boundary conditions. The remaining two surfaces are free of 
applied stress . Numerical work is presented to show the characteristics of 
the solutions and it is worthwhile to note that the stresses are smooth func- 
tions throughout the solid which includes the edges and corners. The solu- 
tions are not general enough to satisfy an arbitrary distribution of applied 
stress. It is recommended that further study be devoted to investigating sums 
of the solutions obtained which may describe additional useful cases. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field , Calif., 94035, May 29, 1967 
124 - 08 - 06 - 01 - 00-21 
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APPENDIX A 

DISCUSSION OF COMPATIBILITY EQUATIONS 


Since the linear compatibility equations govern the small deformations of 
any solid, it is worthwile to see how the relationships between the equations 
can be used without restricting the strains to any particular stress -strain 
law. Reference 6 shows that if all six compatibility equations are satisfied 
on the boundary of a simply connected solid, then satisfying three unrelated 
compatibility equations throughout the interior of the solid is sufficient to 
guarantee continuous displacements. In this section, conditions under which 
three compatibility equations alone are sufficient will be investigated. From 
the following definitions, 


<Pi 


cp 2 


q> 3 


^4 

<P 6 


s 2 e x a 

by 2 4 bx 2 

a r 

'xy 
bx by 



(Al) 

a e x a € z 
az 2 + ax 2 

b^y 

'xz 

dx bz 



(A2) 

fit + V 

bz 2 by 2 

*\ z 
by bz 



(A3) 

2 _ a 

by bz bx 

( _ 

y bx 

by 

+ ay 


(A4) 

2 V . a 



V\ 

(A5) 

bx bz by 

Sx 

ay + 

bz J 


TV + 

^xz 

S V\ 

(A 6) 

bx by bz 

v ax 

ay 

bz J 


the six compatibility equations are then 


<Pj. “ 0 , 

cp 2 = 0 , 

cp 3 = 0 

(AT) 

<P 4 = 0 f 

cp 5 = 0 , 

9 S = 0 

(A8) 


The identities discussed in references 1 to 3 that interrelate equations (Al) 
to (a 6) are 


ckp^ 5cp 5 Sep 
by + bx ~ ^ bz 


(A9) 
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^4 ^ _ 0 5c Pa 

Sz Sx Sy 

acp 5 Sep 6 Scp 3 

Sz by Sx 

or an alternative form is 


S 2 epi 

S ep 2 

S cp 3 

^2 

s cp 4 

8 z 2 

Sy 2 

dx 2 

by Sz 

2 

^2 


2 

0 cpi 

S ep 2 

. . a ^3 

0 9s 

Sz 2 

Sy 2 

8x 2 

bx bz 

2 

. 2 


.2 

s cpi 

C \1 

©■ 

'O 

a cp 3 

d cp 6 

bz 2 

Sy 2 

+ dx 2 

bx by 


(A10) 

(All) 


(A12) 


(A13) 


(A14) 


Since cp , cp 5 , and cp 6 are so closely related to cp^ < 72 , and cp 3 , it seems 
possible that under certain conditions, strain components that satisfy three 
compatibility equations -will automatically satisfy the other three equations. 
To explore this, assume that the strain components satisfy the three 
compatibility equations 


9, = 0 , cp= = 0 , 


9s = 0 


Substituting equations (A15) into equations (A9) to (All) gives 


a<Pi 


Sz ^ ’ by 
Integrating equations (Al 6 ) gives 


Sep, 


= 0 , 


Sep, 

Sx 


= 0 


ep 2 = 


S 2 g x S 2 e y 

s 7xy 

by 2 + bx 2 

bx by 

^2 X 2 

a e x a e 2 

^ 7 X z 

bz 2 bx 2 

bx dz 

^ £y ^ € 2 

a 2 / 

7 yz 

" 


Si= 


(A15) 


(A16) 


(A17) 


Equations (AI 7 ) show that the strains which satisfy equations (A8) may not 
satisfy equations (A7)j "but the residues, at most, are the functions f, g, 
and h. Different forms of the strains can be studied to see what forms can or 
cannot contribute to the residue. It follows from equations (AI 7 ) that some 
forms that could contribute are 
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y 


7 


xy 


fi(x,y) + g x (x,z) 
f2(x,y) + g z (y,z) 
f 3 (x,y) + g 3 (y,z) 

^(x^y) > 7 XZ — ^5(XjZ) j 


"N 


7 yz = fs(y^) 


(Al8) 


In order to find forms that cannot contribute to the residue and, therefore, 
will satisfy equations (A7)> consider those strains that can be represented by 
a power series of the form 


z 


k m n 


G x = / a kmn x y z 


k,m,n=o 


k m n 
b kmn x y z 


k,m,n=o 


G. = 


oo 

z 


k m n 
c kma x y z 


k,m,n=o 




xy 


Z k m n 
a kmn x y z 


k,m,n=o 


xz 


£ Pkmn^z 11 > (A19) 


k,m,n=o 


^yz 


V^J 

Z 


k m n 
?kmn x y z 


k,m,n=o 




where 


OQ 00 00 00 

1 ^ 1 ^ I I 

k,m,n=o k=o m=o n=o 


Equations (Al) to (A3) become 
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k,m,n=o 


CO 


^ ta k ^ m ^^ n (m+l)(m+2) 

k^m^n^o 


+ ^k+2 ,m,n(k+l) (k+2)^ 


- a k+i,m+i,n< k+1 ) ( m+1 > 


cp 2 = 


I 


9 


2 kmn 


t a k,m,n +2 ( n+ -'-) ( n+ 2) + ^.^^^(k+l) (k+2) 


k,m,n=o 


k,m,n=o 

-Pk+i,m,n+x( k+1 ) ( n+1 ) 3 x Vz n 


> (A20) 


CD = 
^3 


z 


9 


3 kmn 


^ ^k,m,n^( n+1 )( n+2 ) + c k,m^,n( m+1 K m+2 ) 


k^m^n=o 


k,m,n=o 

-7k,m+i,n+i< m+l) ( n+l) JxV 1 ^ 11 




wh.ei*e “fcti6 indices h.8/ve t>een changed so "that "tentus of like exponents can he 
grouped. Since the strains given by equations (A19) are assumed to satisfy 
equations (A8) , equations (Al6) are satisfied and can be written as 


my 

kx -1 


^ikmn 

^2kmn 

^ 3 kmn 



(A21) 


Equations (A2l) must be satisfied for all values of x, y, and z so that 

* ^ 


CD 

^-kmn 

= 0 

k = 0 jl j2 y • • • y 

m = 0,1,2, . • . ; 

n “ 1 y2 y • • • 

cp 

^kmn 

= 0 

k = 0 jl j2 y • • • y 

m = 1,2, . . . ; 

n = 0,1,2, . . . > 

cp 

3 kmn 

= 0 

k = 1 y2 y • • • j 

m = 0,1,2, . . . ; 

n = 0,1,2, • • • 

J 


(A22) 


2k 



For equations (A7) to be satisfied, it must also be true that 


9 


1 tano 


= 0 


cp = 0 
2 kon 


cp = 0 
3 omn 


(A23) 


From equations (A20), it is seen that equations (A 23) can be satisfied by 
requiring that, for k,m,n = 0,1,2, . . ., 


a k,m+2,o( m+1 )( m+2 ) 

a k,o,n+2 (n+l) ( n+2) 

b o,m,n-te( n+1 Hn+2) 


+ b k+2,m,o^ k+1 Hk+2) 

+ c k+£,o,n^ k+b ^ k+2 ^ 
+ c o,m+a ,n( m+ - L ) ( m+ 2) 


- a t+i,m+i,o( k+1 Hm+l) 

_ Pk+i,o,n+i( k+ -*-) ( n+b ) 

- ?'o,m+i,n+i( m+1 )( n+1 ) 


0 

° 

0 


(A24) 


Equations (A2^) can be satisfied if the coefficients are combined so that the 
sum is zero or by each coefficient a, b, . . -, 7 being zero. An example 
of a strain that satisfies this latter condition would be if the strain e x 
were of the form 


e x = A cos x sin y sin z 

which has a series expansion where a^. = a ]^ 0 ^ n = 0* 

The conclusion to be drawn from this discussion is that if the strains 
are chosen so that their power series expansions exist and the coefficients 
satisfy equations (A24) , then equations (A8) are sufficient to guarantee that 
all six compatibility equations are satisfied. 

One could have assumed first that the strains satisfied equations (A7) 
and then have looked for forms that would also satisfy equations (A8) . By 
the same procedure, equations analogous to equations (A24) can easily be 
found; but since equations (A12) to (Al4) are of second order, there are six 
equations between the coefficients a, b, . . •, 7- 
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APPENDIX B 


DERIVATION OF EQUATIONS 


Since some manipulation is required to write the compatibility equations 
in the form given by equations (6) and ( 7 ), it seems appropriate to show a 
detailed derivation. 

By Hooke's law the strains can be found in terms of the stress functions. 
Substituting equations (4) into equations (2) gives 


€ x e 


fe,zz + "^3,yy “ / ^ - v ( ^!x,zz + fe,xx _ / Y dy + ^i,yy + ^^xx ~ I % 


dz 


2(1 + v) 


xy 


E 


fe,xy 


(Bl) 

(B2) 


The other strains are symmetric with equations (Bl) and (B2) and can be 
obtained by permutation of the coordinates. The strains can be substituted 
into the first compatibility equation given by 


3 2 € v 3% 3 2 7 


'X + ~ 7 


xy _ 


by bx bx by 


= 0 


(B3) 


which gives 

.2 


a 

ay 2 L 


^2,zz + fe,yy + - ^(v 2 ^!) - J Xdx+v^J Ydy+ j Z dz 

ti ; ZZ + + ^3,yy ~ v(v 2 t2) - f Y dy + \ 


bx 


By definition, 


X dx + Z dz 


V 2 !!/. = ilf. + ilr. + 

Y i,xx Y x,yy y i,zz 


= 0 
(b4) 

(B5) 
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Using equation (B5)> one can write equation (b4) 


as 


dy 2 L 


- v 2 ^ - V 2 ts + (te,3oc + ^2,yy + +3, Z z) 


+ X dx - v Y dy + 




a x " 


-V 2 ti + vv 2 ^ - V 2 ^ 


+ (^l^xx + ^l^yy + te ; zz) + / Y dy - v ( / X dx + / Z dz 


= 0 (b6) 


Since (a^ay 2 ) (^ 2 jXX ) = (a /ax 2 ) (^ 2 ,yy) > e_tc • > these terms can he inter- 
changed to give the form, which is the first of equations (6). The next two 
compatibility equations can be obtained by permutation of the coordinates. 
The fourth equation is 


^ e x d ( ^?yz + ^7 X z + ^7 X y 


dy dz dx \ dx dy 


dz 


= 0 


Substituting the strains gives 

a 2 r 


^2,yzzz + ^s,yyyz 


dy dz 


X dx - v Ui 


-,yzzz + ^3,xxyz 


a y a z f YdJ+ ^ x >yyy z + 


f Z dz^) 
dy dz J /_ 


+ 2 ( 1 + v ) ( y xxyz + ^2 y xxyz + , xxyz ) - 0 


or 


d" 


dy dz L 7 


fc, ZZ + fe,yy " ^l,xx + fc.xx + ^3 


,xx 


,xx 


X dx - v V ti - / Y dy - / Z dz 


= 0 (B7) 


Equation (B5) can be used to reduce equation (B7) to the first of equations (6) 
and the remaining equations can be obtained by permutation of the coordinates. 
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